Scattering bright solitons: quantum versus mean-field behavior 
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We investigate scattering bright solitons off a potential using both analytical and numerical meth- 
ods. Our paper focuses on low kinetic energies for which differences between the mean-field descrip- 
tion via the Gross-Pitaevskii equation (GPE) and the quantum behavior are particularly large. On 
the iV-particle quantum level, adding an additional harmonic confinement leads to a simple signa- 
ture to distinguish quantum superpositions from statistical mixtures. While the non-linear character 
of the GPE does not allow quantum superpositions, the splitting of GPE-solitons takes place only 
partially. When the potential strength is increased, the fraction of the soliton which is transmitted 
or reflected jumps non-continuously. We explain these jumps via energy-conservation and interpret 
them as indications for quantum superpositions on the JV-particle level. On the GPE-level, we also 
investigate the transition from this stepwise behavior to the continuous case. 
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I. INTRODUCTION 

Bright solitons generated from attractively interacting 
ultra-cold atoms have been realized experimentally both 
in quasi one-dimensional (ID) configurations [TJ [5] and 
three dimensions (3D) 

Some of the effects investigated for bright solitons [4] 
could, in principle, also be investigated in classical sys- 
tems. However, bright quantum-matter-wave solitons arc 
mesoscopic quantum objects which are particularly use- 
ful to investigate beyond mean-field effects of quantum 
solitons: While there are cases for which both the classi- 
cal and quantum descriptions agree for particle numbers 
as low as N 3 [5], it has been proposed to use scat- 
tering bright solitons off a potential in order to produce 
non-classical quantum superpositions [5] [7J . 

Other investigations of bright solitons include sta- 
tionary solutions of the (mean field) Gross-Pitaevskii 
equation (GPE) [5J, soliton trains [TO], solitons un- 
der transverse confinement [11) and incoherent matter- 
wave solitons |12j as well as deviation from one dimen- 
sionality [13] and regular and chaotic dynamics in soli- 
ton collisions [13] • Research also covers topics ranging 
from fragmented states [TS] , stabilization and destabiliza- 
tion of second-order solitons against perturbations [T5] 
and quantum reflections [17] over soliton localization 
via disorder [T5J [Tj|] and in time-dependent traps with 
time-dependent scattering length [50], resonant trapping 
through a quantum well [21] and possible applications to 
interferometry [22H24] . 

The focus of the present paper lies on scattering bright 
solitons off a potential in a one-dimensional geometry. 
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On the ./V-particle quantum level this can result in quan- 
tum superpositions; however, the detection of such super- 
positions requires the identification of a clear experimen- 
tal signature. For mean-field solitons both behavior sim- 
ilar to the case of a single particle (see, e.g., Ref. [2"3"ll2"4"] ) 
and considerably deviating behavior has been reported in 
the literature Q By focusing on low kinetic energies, we 
will discuss the parameter regimes for which both types 
of behavior can be expected. 

The paper is organized as follows: Section[Il]introduces 
the models used to investigate scattering bright solitons 
off scattering potentials. The energetically allowed final 
states are discussed in Sec. IHII 

In Sec. |IV[ scattering bright solitons off a scattering 
potential is investigated with additional harmonic con- 
finement as used in experiment [26] . Section |v| explains 
the behavior on the level of the GPE. Section IVII con- 
cludes the paper. 



II. MODELS 

On the level of the ./V-particle quantum mechan- 
ics, the system can be described by the Lieb-Liniger(- 
McGuire) [27] [25J Hamiltonian with additional external 



Scattering mean-field (GPE) solitons off a barrier with behav- 
ior that considerably deviates from a single particle has been 
investigated both in Ref. 1251 and on page 2 of Ref. [7]. 
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with giD < 0. 

The mean-field equation corresponding to the Hamil- 
tonian is the GPE g] 

h 2 

ihd t <p(x,t) =- —d x (p(x, t) + V e xt(x)<p(x,t) 
2m 

+ (N-l)g 1D \^(x,t)\ 2 ip(x,t) , (2) 

where |cp(cc, t) | 2 is the single-particle density (cf. 29). 
While the GPE ^ often is used to describe Bose- Einstein 
condensates (BECs) of finite particle numbers, it can 
strictly speaking only be valid in the limit [29] : 
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There even are cases for which the Gross-Pitaevskii func- 
tional becomes exact in this limit [29] . 

For attractive interactions and without an external po- 
tential [V cxt (a:) = 0], Eq. ^ has exact soliton solutions 
of the form [3] 
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where u is the velocity and xq the initial position. On 
the one hand, Eq. (El describes a single soliton, for which 
normalizing Eq. (El) to one yields (cf. [50] ) 
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On the other hand, Eq. Q can also describe (well sepa- 
rated) parts of a solution. The sum of two such solutions 
which are, e.g., on both sides of a scattering potential 
corresponds to a fraction of the atoms being on one side 
and the rest of the atoms on the other side (note that 
the widths of the two solitons then depends on which 
fraction of the atoms is on each side). This is, however, 
very different from a quantum superposition of one soli- 
ton containing all atoms being simultaneously on both 
sides (cf. [6]). 

Contrary to the localized mean-field solution Q, 
eigensolutions of the iV-particle Schrodinger equation 
corresponding to the Hamiltonian ([I]) with zero external 
potential V^xt have to be translationally invariant (up to 
a phase factor): 
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The eigenenergy E = E (N) + E kin of the eigenfunc- 
tion Q consists of the ground-state energy 
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and the center-of-mass kinetic energy 
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On the one hand, there are cases for which mean-field 
and TV-particle physics can be shown to agree: By as- 
suming the center-of-mass part of the wave function to 
be a delta function, 

tPc(X) = S(X - x - ut) , 

it is possible to show that the single-particle density of 
the Gross-Pitaevskii soliton Q and the many-particle 
soliton ^ coincide in the limit N ^> 1 [31]. Thus, rather 
than simply being the solution of an approximated equa- 
tion, the single-particle density of exact many-particle 
solitons ([6| is well described by the mean- field case Q. 

On the other hand, there are cases for which there is 
no equivalent on the mean-field (GPE) level: Mesoscopic 
quantum superpositions. 



III. MESOSCOPIC QUANTUM 
SUPERPOSITIONS (MQS) VS. 
HARTREE-PRODUCT STATES 

As a starting point to investigate mesoscopic quan- 
tum superpositions (MQS), this section first investigates 
parameter regimes accessible to the A^-particle physics 
when scattering a bright quantum soliton off a poten- 
tial by looking at the energetically allowed final states in 
the iV-particle case (Sec. Ill A) after the soliton has left 
the potential. Section [ill B investigates Hartree-product 
states of the form 
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A. iV-particle quantum mechanics 

It has been suggested that scattering slow bright quan- 
tum solitons of N ps 100 particles off a scattering poten- 
tial generates mesoscopic quantum superpositions [5J [7] . 
Scattering would produce a quantum superposition of all 
particles being either on one side of the scattering poten- 
tial or at the other side: 
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FIG. 1. (Color online) Energetically (dis) allowed states for 
TV = 100 particles with n particles on one side of the scatter- 
ing potential and TV — n particles on other side. The parame- 
ter regions are displayed as a function of n and center-of-mass 
kinetic energy E\& n . (a) Magenta/gray region: not energet- 
ically allowed as soon as all particles have left the scatter- 
ing potential, black and white regions: energetically allowed. 
Only in the white region [which lies above -Bki n — 0.755|i?o| 
for TV = 100; Eq. |j~4| ] are product states with 50:50 occu- 
pation on both side of the barrier energetically allowed, (b) 
Enhanced lower part of the left panel. In the energy regime 
below the dashed line, the only energetically allowed states 
are: all particles on the left or all particles on the right of the 
barrier - and quantum superpositions of both. 



where the Fock-state notation |TV — n, n) denotes TV — n 
particles on the right and n particles on the left of 
the scattering potential. Such states have been called 
"Schrodinger-cat states" or "NOON states" ; further sug- 
gestions how mesoscopic quantum superpositions might 
be obtained can be found, e.g., in Refs. 



MQS can be more general than Eq. ( 11 ); quantum su- 
perposition involving states like 
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or superpositions thereof are also interesting. Further- 
more, quantum superpositions which do not correspond 
to exact 50:50 splitting are also in the focus of research 
(cf. 0). 

Figure [T] depicts the possible energy regimes. For high 
center-of-mass kinetic energy [well in the white area of 
Fig. [l] (a), the parameter regime investigated, e.g., in 
Refs. 
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the soliton is energetically allowed to break into (at least) 
2 parts. The scattering potential can act as a beam split- 
ter on the level of single particles. In this energy regime, 
the mean-field soliton splits [231 121] > which includes 50:50 
splitting. In the high energy regime, the particle radia- 
tion observable at lower energies disappears |43j and the 
system thus stays condensed. Replacing the "^>" by "=" 



white and the black region in Fig. [I] (a): 

, -± -0.755SO : TV = 100 
iiki„ i = _ 75Eq . N 
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For low center-of-mass kinetic energies [below the 
dashed line in Fig. [I] (b)], 



E kin < E (N - 1) - E (N) , 



(15) 



the soliton is energetically forbidden to break into two 
parts and thus all particles are either on one side of the 
scattering potential (n = 0) or on the other (n = TV); 
quantum superpositions of both (111 have been pre- 
dicted [5J. 

In between the threshold given by Eqs. ( 14 ) and 



in Eq. ( 13 1 yields the value for the boundary between the 



(15), contributions of more states are allowed; the ma- 
genta/gray region indicates which |TV — n,n) are ener- 
getically disallowed for each value of the center-of-mass 
kinetic energy. Contributions of states like |TV/2,TV/2) 
automatically imply that this part of the wave function 
has a higher energy than energetically allowed for any pa- 
rameter regime which is labeled magenta/gray in Fig. [T] 
This statement is independent of the TV-particle state 
these particles are in as long as they cannot be found at 
the potential. 

On the TV-particle level governed by the Hamilto- 
nian ([I]), energy is conserved not only in the sense that 
(ip(t)\H\ip(t)) is conserved but also in the sense that all 
higher moments (ip(t)\H v \ilj(t)), v = 2,3,... are time- 
independenlj^ Thus, within the magenta/gray param- 
eter regime depicted in Fig. [T] energy conservation not 
only prevents the final state from being identical to 
|TV/2,TV/2), but |TV/2,TV/2) (and, \N/2-n, N/2+n) with 
increasing n for decreasing initial center-of-mass kinetic 
energy) cannot even be an important contribution to the 
final wave function. 

Just because a certain value for n is energetically al- 
lowed does not necessarily imply that it will occur: at the 
threshold given by Eq. ( p~4j ) , n = TV/2 could, in principle 
occur. However, this would imply that in the final state, 
the fractions of the soliton on both sides of the barrier 
no longer move. It is thus more likely to occur for even 
higher initial center-of-mass kinetic energies. 



B. Mean-field approach via Hartree-product states 

Often, mean-field theories are introduced [4] for 
bosonic TV-particle quantum systems by starting with 



2 As the time-evolution operator U(t,0) = exp(—iHt/h) (with 
|V(*)> = U(t,0)\tp(0))) commutes with H" {v = 1,2,3,...), we 
have: (iP{t)\H v \%p{t)) = (%p(0)\H v \^{Q)) for v = 2, 3, 4, . . .. This 
excludes contributions from eigenfunctions for which the eigenen- 
ergy is not negative enough to contribute to the final state. Alter- 
natively, H" could be replaced by [H — Eq{N)] li before repeating 
this analysis. 



4 



Hartree-product states ( 10 ), which are subsequently used 
to derive Gross-Pitaevskii equations like Eq. ^ [4 . This 
does, however, not necessarily imply that GPE is equiva- 
lent to the Hartree-product states; rather than interpret- 
ing the GPE-solution <fi(x, t) automatically as being part 
of a Hartree-wave function, a more general approach is 
to regard \<fi(x,t)\ 2 as the "single-particle density" [25] . 

Hartree-product states for which (p(x,t) is zero on one 
side of the potential always exist. More interesting to 
compare with Sec. |III A| are wave functions for which 
<p(x, i) is non-zero on both sides of the scattering po- 
tential. In this case, both the Fock-state \N/2, N/2) and 
\N/2 + n, N/2 — n) (with small n) are thus involved in the 
many-particle Hartree-product wave function: All Fock 
states contribute to the wave-function [as long as the 
single-particle wave function is non-zero at both sides, 
cf. Eq. (A9)]. Note that contrary to the case of Hartree- 



product states, contributions of such states to the to- 
tal wave function can be avoided in the case of a MQS. 
Strictly speaking, for low center-of-mass kinetic energies, 
those states are not accessible on the Af-particle level (as 
discussed in the previous section). This seems to agree 
with the rather stepwise behavior of scattering mean-field 
solitons reported in Refs. [3 US]. This paper investigates 
the transition from this stepwise behavior to the continu- 
ous case similar to single-particle physics in more detail. 
We will show that Fig. [T] provides the energy scale on 
which the non-splitting GPE-soliton becomes a splitting 
soliton, depending on its initial center-of-mass energy. 
The technical details are presented in Appendix |A") 



IV. SCATTERING DYNAMICS IN THE 
PRESENCE OF HARMONIC CONFINEMENT 



A recent (quasi-) ID experiment 26J so far com- 
bines scattering solitons off a potential in ID with addi- 
tional harmonic confinement. We model this situation in 
Sec. IV A So far, it does not yet realize the regime of both 
low kinetic energies and temperatures necessary to pro- 
duce the quantum superpositions suggested in Refs. [S][7]. 



For the A^-particle quantum case (Sec. IV B), the idea 
is to start with the many-particle ground state in the 
harmonic trap (cf. [44]), the center of the trap can then 
(quasi-)instantaneously be shifted, followed by switching 
on the scattering potential in the middle of the trap. 
For the mean-field description via the GPE (Sec. IV C) 



the same situation is repeated for a single soliton with- 
out initial center-of-mass kinetic energy. While the 
two Sees. |IVB| and |IV C| cover different energy regimes, 
Sec. IV D demonstrates what effects can happen, in prin- 
ciple, on the Af-particle level for the parameter-regime for 
which the GPE displays the peculiar behavior discussed 
in Sec. HVCl 



A. Harmonic confinement and scattering potential 

We assume the scattering potential to be narrow 
enough^] for it to be approximated by a delta function: 



V cxt {x) = -Wi 2 + v 5(x) ; 



(16) 



where "narrow" refers to the potential being narrower 
than both the soliton width (cf. [HI US]) and the oscilla- 
tor length: 



Agpe 



(17) 



where the index GPE indicates that this is a relevant 
length scale for the GPE. 



B. iV-particle quantum physics: effective potential 
approach 



The effective potential approach of Ref. (6[ [45] is valid 
for low center-of-mass kinetic energies^] the system can 
be described by an effective Schrodinger equation for the 
center-of-mass coordinate X. This effective Schrodinger 
equation reads: 
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+ V cff (A)Vc(A,i) 
where the effective potential is given by j6[ 05] 

V eS {X) = J d N x\i/j Ntk (x)\ 2 V(x)6 (x-^E^Li^) 
U 



cosh 2 (A/^) ' 

where the last line uses the results of Ref. 
parameters are: the strength 



U 



Nv m\g lu \(N - 1) 

4 h 2 



(19) 
and the 

(20) 



which is the product of the particle number, the single- 
particle potential strength v [Eq. ( 16 1] and the soliton 
amplitude [Eqs. Q and and the width 



(21) 



m\g 1T) \(N-l) 



3 Wide potentials were discussed in the TV-particle quantum case 
in Refs. OE]. 

4 The validity of the effective-potential approach \E\ 1451 to de- 
scribe time-dependent scattering has been proved rigorously by 
calculating strict bounds on the transmission and reflection am- 
plitudes [46]. 
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given by the soliton width JUI, cf. Eq. 

The effective potential ( 19 1 thus has the form of the 
soliton; the ratio of the width I to the center-of-mass 



oscillator length [cf. Eqs. (17 1 and (18)] 
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^GPE 
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influences the physics. 

As soon as the center-of-mass density near the effective 
potential approaches zero, MQS states emerge (cf. [5]). 
For center-of-mass energies which are an order of magni- 
tude higher, numerics still predicts quantum superposi- 
tions far from product states [7J. 

One advantage of the additional harmonic trapping po- 
tential is that excitations due to an opening of the trap 
which prepares the initial state jUJ can be discarded. 
Another advantage is discussed in this section: simply 
producing potentially interesting MQS and showing that 
in an experiment all particles are either on one side of 
the potential or at the other would not be enough to 
confirm the existence of MQS; additional experimental 
signatures are necessary. One possibility is to look at the 
interference in the center-of-mass density after removing 
the barrier [5J. 



For not-too-broad effective potentials (19), there is a 
simpler approach shown in Fig. [2] After the second col- 
lision of both parts of the wave function in the presence 
of the barrier, this leads to a probability close to one of 
all particles being at the side opposite to the initial con- 
dition [Fig. [2] (b)], which considerably differs from the 
case of a statistical mixture [Fig. [2] (a)]. Losing a single 
particle after the creation of the quantum superposition 
would turn the MQS (11) into such a mixture. While 



the precise influence of decoherence will depend on ex- 
perimental details, we model decoherence by replacing 
the quantum superposition of all atoms being either on 
one side of the scattering potential or on the other by a 
statistical mixture near t ss T/2. 

In order to understand this behavior, let us first assume 
the scattering potential is narrow enough to be approxi- 
mated by a delta function: 



VciTpO = —m(x + x s ). 

m 



(24) 



where we included a small shift to the left. 

In order to derive the leading order behavior of our 
center-of-mass wave packet, start with the fact that with- 
out the scattering potential, after some time the wave 
packet would be in the middle of the potential with all 
the initial potential energy being transformed into ki- 
netic energy. For a plane wave with wavevector k, the 
transmission coefficient is given by 49J 
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ik 



(25) 




X/ \ oi 



FIG. 2. (Color online) Center-of-mass density of a ID bright 
quantum-matter-wave soliton in a two-dimensional projection 
as a function of both time t (in units of the oscillation period 
T) and elongation x (in units of the oscillator length A osc ). 
There is additional harmonic confinement (oc x 2 )\ the soli- 
ton is modeled within the effective potential approach. Ini- 
tially, the wave function is centered around x = — 10A OSC . 
At t — 0.25T the soliton scatters off a narrow effective po- 
tential modeled by a delta potential and the center-of-mass 
wave function splits into two parts, leading to a mesoscopic 
quantum superposition of all particles either being on the left 
or right of the scattering potential at t — 0.5T. (a) If the 
quantum superposition becomes a statistical mixture due to 
decoherence at t ~ 0.5T, there is a 50% probability to find the 
particles on either side of the scattering potential at t = T. 
(b) For a quantum superposition, all particles end on the 
right side with p T = 98.5% probability at t — T. 



and reflection coefficient by [49] 
R ^ 



ik — Q 



(26) 



For 50:50 splitting, they must have the same modulus, 
which thus defines the strength of the potential: 



(27) 



In our case, we do not have plane waves but wave packets 
centered around k — fl. The time scales for the reflection 
of both parts of the wave packet in the harmonic potential 
are the same. 

Due to the harmonic confinement, in the further time 
evolution a second scattering takes place where it has to 
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FIG. 3. (Color online) (a) For the same situation as in Fig. [2] 
the probability to find all particles on the right side is plot- 
ted as a function of the width of the effective potential \19) . 
The cases of the quantum superposition (solid line) and the 
statistical mixture (dashed line) are clearly distinguishable 
for repeated measurements. Points: Averaging over an ex- 
perimentally realistic narrow Gaussian distribution with 
width 5 centered around N = 100 and truncated at 90 and 
110. (b) Probability p r to find all particles on the right of 
the scattering potential as a function of shift x 3 . From top 
to bottom: £/X OBC = 0,0.125,0.5,1.0. (c) Difference of nu- 
meric solution and the computer-algebra based analytic so- 
lution extending ( 29 1 to general potentials ( 19 1 for the same 



parameters as in the previous panel, (d) Probability to find 
all particles on the right hand side as a function of the shift 
for £ — A osc . Compared to Eq. ( 29 1 , the amplitude is lower 
than one (caused by the wider potential). Furthermore, the 
fact that the shift of the potential leads to an interference of 
only parts of the wave packets [discussed below Eq. ( |29[ l] is 
clearly visible. 



be considered that the transmitted part had to cover an 
additional distance of Ax s : 



u(x) 



(28) 



The transmission coefficient after two reflections is there- 
fore given by 



T 



4iUx s 



(i-ir (i-ir 

- [1 + cos(4ftx s )] . 



(29) 



Thus, if the scattering potential can be approximated by 
a delta function, we can indeed expect a probability close 
to one that all particles are on the right side of the scat- 
tering potential if they initially were on the other side. 



A more complete analysis would have to include wave 
packets rather than plane waves. This will effectively 
damp the oscillation amplitude as a function of the dis- 
placement x s as can be seen in Fig.|3](d). In this panel, 
the probability to find all particles on the right does not 
reach 1 because a broader effective potential (cf. [6]) was 
used. 

In Fig. [3] we display the probability to find all particles 
on one side of the potential after one oscillator period. 
The quantum mechanics of pure states can clearly be dis- 
tinguished from statistical mixtures over a wide param- 
eter regime. The fact that the above scheme depends on 
the position of the scattering potential offers a potential 
application to interferometry by identifying small poten- 
tial gradients along the center of the harmonic trap (for 
applications of ultra-cold atoms in interferometric experi- 
ments discussed in the literature see, e.g., [23, 24, 50, 51 ). 

While the difference between pure quantum dynam- 
ics and statistical mixture is particularly large for a very 
narrow scattering potential, it is still clearly visible for 
broader effective potentials [Fig. [3] . The two values are 
distinguishable as soon as the error of the means are small 
enough (which scale as 1/ y/N Icp , where N lep is the num- 
ber of repetitions of the experiment) . The experimentally 
realizable parameters for a soliton of N — 100 particles 
discussed in Ref. [5] correspond to I = 1.5A osc - where 
the difference between statistical mixture and MQS is 
still clearly visible. 



C. Mean- field approach via the GPE 



Modeling the same situation as in Sec. IV B on the level 
of the GPE leads to a different energy regime because of 
the mean-field limit In this limit, the effective po- 
tential regime below the dashed line in Fig. [I] (b) could 
only be covered for vanishing ratio of center-of-mass ki- 
netic energy to ground-state energy. However, the GPE 
can cover the regime for which this ratio is finite. 

Based on the reasoning of Sec. IIIB (cf. Sec. A 3 1, it 



would not be surprising to find a stepwise behavior of the 
reflection (or transmission) coefficient which jumps from 
to 1. 

Figure [4] shows a different behavior: while there are 
indeed jumps in the reflection coefficient if the strength 
of the scattering potential is increased, this jump lies be- 
low 1 for many parameters and gets smaller if the kinetic 
energy approaches the threshold (14). When repeating 



the calculation without harmonic confinement, the qual- 
itative behavior in all three cases is the same. 

Furthermore, changing from energies for which the 
product state corresponding to 50:50-splitting cannot ex- 
ist on both sides of the scattering potential [Fig. E] (a)] 
to energies for which it can (just) exist [Fig. [4] (bj\ pri- 
marily reduces the size of the gap (for Fig. [4] by a factor 
of 2) before it eventually vanishes for even higher kinetic 
energies. In order to investigate this in more detail, the 
next section focuses on a more detailed analysis without 
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V^kin^GPE) 

FIG. 4. (Color online) Reflection of a GPE-soliton in a ID 
harmonic trap from a narrow Gaussian barrier [used to model 
the delta-potential in Eq. ( 19 1] as a function of strength of 
the scattering potential. The initial center-of-mass kinetic 
energy increases from top to bottom, (a) Emu = 0.7|£Jo|: 
The reflection R remains below 0.22 before jumping to values 
above 0.77 near v / (E^Xgpb) = 2.65126. (b) E kin = 0.8|£ |: 
The reflection R remains below 0.36 before jumping to values 
above 0.64. (c) Skin = 1.0|£o|: The reflection 7? is continuous. 



where Cj are the boson creation and annihilation oper- 
ators on site j, fij = c\c- are the number operators, J 
is the hopping matrix element and U the on-site inter- 
action energy. For a small Hilbert space, such a model 
can be solved via exact diagonalization (see, e.g, and 
references therein); for a larger Hilbert space, imaginary 
time evolution is a much better choice to determine the 
ground state, that is our initial condition (cf. [53]). We 
use this to find our initial condition and subsequently 
numerically solve the full time-dependent Schrodinger 



equation corresponding to the Hamiltonian ( 30 1 via the 
Shampine-Gordon routine |54j . 

While the limit ^ is not accessible on the ./V-particle 
level, N = 4 still allows us to use the full Hamiltonian and 
get physical insight into what happens on the iV-particle 
quantum level. For future quantitative comparison with 
experiments, advanced approximate methods on the N- 
particle level as used in Refs. [3 (53j will be useful. 

Figure [5] shows the difference between the regime of 
perfect MQS like ( [TT|) in the regime of low center-of-mass 
kinetic energy [Fig. IsJ (a), cf. Sec. |IVB] compared with 
the high kinetic energy regime [FigJ5j (c)] and the regime 
of medium kinetic energies [Fig. pKb)]. In the high- 
energy regime, the final state corresponds to a product 
state for which the probability distribution of particle 
numbers always has a single peak. Bimodal distribu- 
tions correspond to quantum superpositions of (for small 
particle numbers) primarily two states, |m, N — n\) and 
\ri2, N—n2) for which the n\ and n 2 differ by more than 1. 



the harmonic confinement. 



V. BOUNDS AND QUANTITATIVE 
PREDICTIONS FOR THE GPE-SOLUTIONS 



D. iV-particle quantum physics: beyond the 
effective potential approach 



One approach to include more particles is to discretize 
the iV-particle Schrodinger equation corresponding to the 
Hamiltonian ([!]) , which is a delicate matter for attractive 



Without the harmonic confinement, the qualitative be- 



havior is the same as in Sec. IV C in addition, as soon as 



the probability to find particles on the scattering poten- 
tial is small, the exact eigensolutions of the Lieb-Liniger 
equations are sufficient to expand the wave function. 

Including center-of-mass kinetic energy and using 
again the notation that n particles are on one side and 
JV-non the other side of the scattering potential, we 






FIG. 5. (Color online) Probability to find n particles on 
the left side of the scattering potential as a function of time 
r = Jt/h for a total particle num ber of N = 4. Initially, 
the system [modeled by Eq. (30 1 on 51 lattice sites] is in 
the iV-particle ground state (determined via imaginary time- 
evolution) of a one-dimensional harmonic oscillator shifted 
by 12 lattice sites before removing the shift and switching on 
the scattering potential. With strong attractive interactions 
[NU/J — —4, panel (a)], the probability to find 4 particles 
on either side of the scattering potential lies above 98% for 
r « 50. For weaker interactions [NU/J — —3, panel (b)] this 
value has already dropped to less than 78%. Without interac- 
tion [NU/J = 0, panel (c)], the value has dropped to 12.5%, 
the value predicted for product states [see Eq. ( A9 1] . The 
values for the scattering potential (vq/J = {0.27,0.43,0.82}) 
were chosen via bisection; A — 0.001J. 



can use energy conservation^] 

£ ki „ + Eo(N) = £ km ,2 + E (n) + E (N - n). (31) 



5 We again assume that the wave function vanishes at the po- 
tential. This allows us to use the exact eigenfunctions of the 



While the left-hand side describes the initial condition, 
the final situation does not necessarily have to consist 
of exactly two solitons as the soliton is energetically al- 
lowed to break into more than two parts. This does 
not influence the exact bounds we derive in this sec- 
tion but it might affect quantitative predictions derived 
from Eq. (31 1. While this section does not explicitly in- 
volve Hartree-product states, Appendix [X] shows that a 
Hartree product state leads to replacing the variable n 
in Eq. (31) by its mean n. After this replacement, the 



derivations of this section remain valid even for Hartree- 
product states. 

Dividing Eq. (31) by \Eq(N)\, and taking the mean- 
field limit ([3]) yields: 



-E<kin — Bkin,2 

\Eo(N)\ 



+1 - R 3 - (1 - R) 3 

3R{1 - R) (32) 



where 



N 



(33) 



is the fraction of the particles which are reflected. Note 
that this argument assumes that the soliton breaks into 
exactly two parts - for more than two parts the final 
center-of-mass kinetic energy -Ekin,2 would be even lower. 

As the final kinetic energy cannot be negative, Eq. ( 32 ) 
implies 



^kin 

\E (N)\ 



> 3(1 -R)R 



(34) 



Thus, for small kinetic energies not all values for R are 
allowed. To proceed, we put all R into the sets {R < 0.5} 
and {R > 0.5} (or in both in case the value R = 0.5 
occurs). Defining 



and 



R« 



Ru 



m.in{R > 0.5} 



max{i? < 0.5} 



(35) 



(36) 



implies that for very low center-of-mass kinetic energy, 
R « 1 or R « are possible - contrary to the quantum 
case which (depending on the strength of the scattering 
potential) allows all values for R in the range < R < 1. 

To test this statement, the magenta/gray region in 
Fig. [6] (a) shows the energetically forbidden parameter 
regime for R as described by Eq. p4| . Numerically cal- 
culating the values for i? max and i? m i„ yields that they 
indeed lie outside this regime. The jumps of the reflec- 
tion reported in Refs. [7j[25] are also visible in this figure; 
we can quantify that it takes place for center-of-mass ki- 
netic energies below one fourth of the modulus of the 
ground-state energy. 

In order to test the prediction of Eq. (31), the center- 



Lieb-Lingcr model, cf. Eq. (120) of Ref. [30] 



of-mass kinetic energy of the solitons after scattering off 
the potential has to be evaluated numerically. To do this, 
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VI. CONCLUSION 




B kir/l £ ol 



0.5 



(b) 



0.5 



1.5 



v /v c 



FIG. 6. Scattering a mean-field soliton off a delta potential 
in a ID situation without additional harmonic confinement, 
(a) The reflection 7? m in (solid line) and i? m ax (dotted line) as 
defined in Eqs. ( 35 1 and ( 36 1 as a function of initial center-of- 



mass kinetic energy. As predicted by Eq. (34 1, both curves lie 
outside the magenta/gray region for energetic reasons, (b) 
Scattering of a GPE soliton with initial center-of-mass ki- 
netic energy i?kin = 0.5|i?o| as a function of the strength 
of the delta-scattering potential, normalized such that with- 
out interaction, there would be 50:50 splitting [cf. Eq. \27\ ]. 
From top to bottom: Dotted red/black line corresponds to the 
the fraction p so iit on of the final GPE-solution which is in soli- 
tons (the rest being particle-radiation). Solid magenta/gray 
line: Total kinetic energy after the scattering as predicted 
by Eq. ( |32[ | using the numerically calculated values for R. 
Dashed blue/black line: numerically calculated total kinetic 
energy. The difference between the two lower lines is only 
large if p so uton < 1. 



We investigate scattering bright solitons generated 
from potentials both on the mean-field level and on the 
TV-particle quantum level in ID. Adding an additional 
harmonic confinement leads to interesting effects both on 
the ./V-particle quantum level and on the GPE mean-field 
level: 

Firstly, adding a harmonic confinement to the cre- 
ation of quantum superpositions of slow bright quantum 
matter-wave solitons provides a possibility to distinguish 
quantum superpositions from statistical mixtures: after 
scattering off the potential twice, the probability that all 
particles are on the side of the potential opposite to the 
initial condition clearly distinguishes the two cases. As 
on the single-particle level, changing the strength of the 
potential leads to a continuously varying reflection (or 
transmission) coefficient. 

Secondly, for the reflection behavior of Gross- 
Pitaevskii solitons, we have derived analytic bounds on 
the size of the jump of the reflection coefficient and de- 
rived the energy scale on which the step vanishes contin- 
uously for increasing center-of-mass kinetic energy. This 
bridges the two types of behavior previously reported 
in the literature: bright solitons on the Gross-Pitaevskii 
level have been reported to split for larger energies (cf. 
[231 RM] ) and display a stepwise reflection-behavior for 
lower center-of-mass kinetic energies [25] for which the 
GPE-soliton hardly splits. 

While the nonlinear character of the GPE does not al- 
low GPE-solutions which arc quantum superpositions, we 
conjecture that the jumps in the transmission behavior 
of the GPE-reflection coefficient indicates the formation 
of interesting quantum superpositions at the N-particle 
quantum level. 

Note added: the jumps on the GPE-level were re- 
cently also investigated in Ref. [56 . 
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Appendix A: Re-deriving the bounds of Sec. [V] for 
Hartree-product states in the limit N — > oo 



This section deals with Hartree-product states ( 10 1 for 



which the single-particle wave function <p(x,t*) is zero 
near the scattering potential and non-zero on both sides 
of the scattering potential at a certain point of time t* . 
Using properties of the eigenfunctions of the Lieb-Liniger 
model [5U] (Sec. A 1 ), we extend the calculations of Sec.|v| 
to Hartree-product states. While this could, in princi- 
ple, lead to different results, in the limit N — > oo we 

(Sec 



recover the same bounds as before 
by assuming that the higher energy states not accessi- 
ble on the iV-particle quantum level are also inaccessible 
to Hartree-product states (cf. the footnote on page [3]), 
Hartree-product states would predict a different reflec- 



tion behavior than observed at the GPE-level (Sec. A 3 1 



1. Justification why the eigenfunctions of the 
Lieb-Liniger Model can be used 

We take the scattering potential to have the form 

V scatt (x) = for |x| >C, (Al) 

which is fulfilled exactly for the delta function used in 
this paper and approximately for Gaussian or 1/ cosh 2 (a; ) 
potentials. 

The single-particle wave function defining the Hartree- 
product states ( 10 1 is assumed to fulfill: 

= J \ip(x,t*)\ 2 dx . (A2) 

This requirement was fulfilled in the numerics of Sec. |Vj 
in the following we use that 



p = / \<p(x,f)\ 2 dx 



implies 



\<p(x, t*)\ 2 dx = 1 — p . 



(A3) 



(A4) 



Equation ( A2 ) implies that we can use the eigenfunc- 



tions of the Lieb-Liniger model to expand the wave- 
function. If n particles are on one side, they can form 
an n-particle soliton with energy: 

(E) = (£ kin ) + E (n) (A5) 

where Eq is given by Eq. pL Alternatively, it could 
form several solitons or even free particles. However, any 
such combination has a higher energy Eq as Eo(n + n) < 
E (n) + E Q {n). 

Thus, if n particles are on one side of the scattering 
potential and N — n on the other (which we denote by 
\n,N — n}) and using the fact that the kinetic energy is 
non-negative, we have: 



2. Bounds on the energy in the limit N — > oo 



Using the notation introduced in Eq. ( 11 ), we can write 



the total Hartree-product wave function at time t* as: 

N 



(A7) 



71=0 



which yields: 



JV 



A2| . Only (ip(x,t*)\H\ip(x,t*)) = ^ \c n \ 2 (n, N - n\H\n, N - n) 



(A8) 

For Hartree-product states, the \c n \ 2 are given by the 
Binomial distribution [57J 3.2.2]: 



N 



p n (l-p) 



N-n 



(A9) 



where p is given by Eq. ( A3 ) . This distribution is strongly 
peaked around 



n = Np 

with root-mean-square deviation 

a = ^Np(l -p) 



(A10) 



(AH) 



Using Eq. ( A6 ) we thus have 



(il>{x,t*)\H\%l){x,t*)) ~ (n,N -n\H\n,N -n) 

>E a (n)+E (N-n); (A12) 

the larger N, the better the approximation in the first 
line of Eq. (A12) becomes. By replacing Eq. (33) with 



R 



N 



(A13) 



(which is the same as p) , the equations and bound derived 
in Sec. |V| are also valid for Hartree-product states. For 
the relation of the initial center-of-mass energy and the 
reflection coefficient we thus reproduce the equation ( 34 1 : 



-Ekii 



\Eo(N)\ 



> 3(1 - R)R 



(A14) 



This leads again to the statement that the magenta/gray 
region of Fig. [6] is forbidden energetically. 



3. Possible differences between Hartree-product 
states and GPE: energy fluctuations 

The above calculation did, however, use the limit N - 



oo for which the distribution (A9) is very narrow. For 



finite N and low initial center-of-mass kinetic energy [cf. 
Eq. (A14)], energetically disallowed eigenfunctions like 



(n, N - n\H\n, N - n) > E a (n) + E (N - n) (A6) those corresponding to \N/2, N/2) would always be an 
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important contribution to the sum ( A7). The only way to 
prevent this (cf. footnote on page [3) is to predict that all 
particles have to be on one side of the scattering potential 
(i.e., R = or R = 1). 

Hartrce-product states for low particle numbers would 
thus behave differently from the GPE, which is strictly 
speaking only valid in the limit pi). Note that 
the contribution of energetically disallowed states like 
\N/2, N/2) vanishes only in the limit N —> oo for which 
limN^^ \c N / 2 / c n\ = (it N/2 ^ n), although taking this 
limit first and then apply energy arguments might be a 



physically relevant approach. 

Thus, by extending this energy argument to higher 
particle numbers one could construct an example for 
which a Hartree-product state disagrees in its "digital" 
prediction for the reflection coefficient with the GPE. 
This would, however, not confute the GPE. The fact that 
Hartree-product states are often used to derive the GPE 
(cf. [3]) does not necessarily imply that both always are 
equivalent; in general, the solutions of the GPE can al- 
ways be interpreted as |v?(x,£)| 2 being the single-particle 
density [55] rather than automatically being part of the 
Hartree-product state (tTol). 
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